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Abstract 

This paper presents the development of a generic 
component level model of a turbofan engine 
simulation with a digital controller, in an 
advanced graphical simulation environment. The 
goal of this effort is to develop and demonstrate 
a flexible simulation platform for future research 
in propulsion system control and diagnostic 
technology. A previously validated FORTRAN- 
based model of a modem, high-perfomiance, 
military-type turbofan engine is being used 
to validate the platform development. The 
implementation process required the 
development of various innovative procedures, 
which are discussed in the paper. Open-loop and 
closed-loop comparisons are made between the 
two simulations. Future enhancements that are 
to be made to the modular engine simulation are 
summarized. 


Nomenclature 


al4 : 
al6 : 
a8 : 
a9 : 

AB: 
bprl4g : 
CAD: 

CLM: 

HPC: 

HPT: 

1GV: 

LPT: 

MAPSS: 

p2: 

P 27: 


Forward VAB1 area actuator (in 2 ) 
Aft VAB1 area actuator (in 2 ) 
Nozzle throat area actuator (in 2 ) 
Nozzle exit area actuator (in 2 ) 
Afterburner 

Forward duct bypass ratio 
Controller and Actuator 
Dynamics 

Component Level Model 
High-Pressure Compressor 
High-Pressure Turbine 
Inlet Guide Vane 
Low-Pressure Turbine 
Modular Aero-Propulsion 
System Simulation 
Fan inlet pressure (psia) 

HPC pressure (psia) 


PL 

Proportional-plus-integral 

PLA: 

Power Lever Angle (degree) 

pnrm2 : 

Percent corrected fan spool 
speed 

ps3: 

HPC exit static pressure (psia) 

psl5 : 

Bypass duct static pressure at 
mixer (psia) 

ps56 : 

LPT exit static pressure (psia) 

sm2 : 

Fan Stall Margin 

stp2 : 

Fan variable inlet stator vane 
actuator 

stp27 : 

HPC stator vane actuator 

stp27d : 

Booster tip stator vane 
actuator 

t2: 

Fan inlet temperature (°R) 

t27: 

HPC temperature (°R) 

t3: 

HPC exit temperature (°R) 

t56: 

LPT exit temperature (°R) 

t5b: 

LPT blade temperature (°R) 

tmpc\ 

Average core metal temperature 
aft of burner. 

VAB1: 

Variable Area Bypass Injector. 

wl4 : 

Forward duct flow rate (slugs/s) 

w2 : 

Fan exit flow rate (slugs/s) 

wf36\ 

Burner fuel flow actuator (pph) 

wf6: 

Afterburner fuel flow actuator (pph) 

xnh : 

High-pressure spool speed 
(RPM) 

xnl : 

Low-pressure spool speed 
(RPM) 


Introduction 

In recent years, interest in intelligent engine 
technology with advanced active component 
control has increased. This interest is being 
pursued in areas such as life extending control, 
which reduces the effect of engine component 
damage and extends engine life while 
maintaining optimal component performance 
[1]; and model-based control, which integrates a 
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real-time engine model into the control 
architecture [2], The computer engine models 
used in these and other types of control studies 
are component level models (CLM), a model 
that links individual component models of state 
space and non-linear algebraic equations, written 
in a computer language such as FORTRAN. 
The difficulty faced with performing control 
studies on FORTRAN -based models is that the 
FORTRAN platform is not supported with 
control design and analysis tools, and there is no 
means to do real-time control implementation. 
It is desirable to have a simulation environment 
that is straightforward, has modular graphical 
components and allows easy access to health, 
control, and engine parameters through a 
graphical user interface. Such a tool should also 
provide the ability to convert a control design 
into real-time code. Because of the current 
technology in advanced modeling software, such 
as MATLAB r [3] and Simulink® [4], such an 
engine model, represented in a graphical 
simulation environment, has the capability to 
become an extremely powerful tool in control 
and diagnostic system development. 

A generic turbofan engine simulation is being 
implemented as a modular aero-propulsion 
system simulation (MAPSS) using Simulink®, a 
computer aided control design and simulation 
package that allows graphical representation of 
dynamical systems in a block diagram form. 
The goal of this effort is to develop and 
demonstrate a flexible turbofan engine 
simulation platform. This platform will allow 
easy access to health, engine, and control 
parameters, as well as a quick way to test control 
and diagnostic systems. 

This paper focuses on the development of the 
MAPSS environment and its validation using a 
FORTRAN engine simulation. The engine 
model to be simulated is described first. Next, 
the paper discusses the approach and validation 
methods used for implementation, including 
example comparisons between the two 
simulations. Finally, the paper concludes with a 
summary of the validation results and future 
work to be done. 


Engine Model Description 

A previously validated FORTRAN-based model 
of a modem high-pressure ratio, dual-spool, low 
bypass, military-type engine, with a digital 
controller, is being implemented in MAPSS. 
This engine model (seen in Figure 1) is a 
nonlinear, low bandwidth, transient, 
perfomiance model similar to the type used in 
various areas of intelligent engine control 
research such as model-based control [2] and 
nonlinear performance seeking control [5], 

The FORTRAN simulation has an updating rate 
of 50 Hz, which is the rate needed by its 
controller. The engine model is represented as a 
CLM, along with accompanying maps that 
model component characteristics [5]. The 
components of this engine model consist of a 
single stage high-pressure ratio fan with variable 
inlet stator vanes, booster with independent hub 
and tip stator vanes, high-pressure mixed flow 
compressor, double-annular combustor, highl- 
and low-pressure turbines, afterburner, and 
nozzle components. Also included are forward 
blocker doors and an aft variable area bypass 
injector (VAB1), giving the engine variable 
cycle capability [5], The CLM uses low- 
pressure and high-pressure rotor speeds, as well 
as average hot section temperature (i.e. metal 
temperature) as state variables. The integration 
of these state variable occurs within a while 
loop in order to achieve mass, energy, and 
momentum balance within a certain tolerance. 

Along with a diagram of the implemented 
engine model, Figure 1 also lists the actuators 
and sensor variables used by the controller. The 
actuators, listed above the engine diagram, are 
for the fan variable inlet stator vanes ( stp2 ), 
forward blocker door area positioning (a 14), 
booster tip stator vanes ( stp27d ), high-pressure 
compressor and booster hub stators vanes 
(, stp27 ), burner fuel flow (yvf36), aft VABI area 
positioning (a 16), afterburner fuel flow (wf6), 
and nozzle throat and exit area positioning ( a8 
and a9, respectively). There are many sensors 
placed throughout the engine, as shown in 
Figure 1 below engine diagram. The sensors 
used by the controller in this model measure 
fan inlet temperature (t2) and pressure (p2). 
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low-pressure spool speed (xnl), high-pressure 
compressor inlet temperature (t27) and pressure 
(p27), high-pressure spool speed (xnh), high- 
pressure compressor exit temperature ( t3 ) and 
static pressure (ps3 ), low-pressure turbine blade 
temperature (t5b), bypass duct static pressure at 
the mixing plane ( psl5 ), low-pressure turbine 
exit temperature (t56) and static pressure (ps56 ). 

This FORTRAN simulation also has a digital 
controller model which runs at a fixed frame rate 
of 50 Hz for sampling sensor inputs and issuing 
actuator commands. This controller is 
represented in Figure 2. Inputs used by the 
controller are power lever angle (PLA), Mach 
number, altitude, and sensor outputs from the 
CLM. The controller uses an open loop 
scheduler to control the stator vane actuators. 
Proportional-plus-integral (PI) control is used to 
regulate actuator demands for the burner fuel 
flow ( wf36 ), nozzle exit area (a8), and aft VAB1 
exit area ( al6 ) through five modes: high, 
medium, and low PLA power settings, fan stall 
margin ( sm2 ) and overspeed mode. Each mode 
has its own set of gains for PI control. The PLA 
is also used in control logic to position al4. The 
actuator demands set the corresponding actuator 
positions used by the CLM (see Figure 3). The 
CLM then feeds the controller calculated and 
measured sensor values in order to re-compute 
actuator demands. 

The goal of this paper is to explain the process 
used to implement these functions and features 
in the MAPSS environment, and thereby obtain 
a new, powerful development tool for control 
and diagnostic systems. To explain this process, 
the next section gives a description of the 
MAPSS environment, followed by a discussion 
on the implementation approach and validation 
procedures that were employed to develop 
MAPSS. 


MAPSS Description 

The Modular Aero-Propulsion System 
Simulation (MAPSS) is a nonlinear, non-real 
time system composed of the “Controller and 
Actuator Dynamics” (CAD) and “Component 
Level Model” (CLM) modules. MAPSS is 
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capable of aiding in the same areas of research 
as the FORTRAN model with increased 
flexibility for control systems development, as 
discussed before. 

The controller in the CAD module emulates the 
functionality of a digital controller, which has a 
typical update rate of 50 Hz. The CLM module 
simulates the dynamics of the engine 
components and uses an update rate of 2500 Hz, 
which is needed to iterate for solutions for mass- 
energy balancing among system components. 
The actuators in the CAD module use the same 
sampling rate as the CLM. The method used to 
obtain the mass-energy balance will be 
discussed in the next section. 

The graphical user interface (GUI) 
accompanying MAPSS is designed to simplify 
the user’s interaction with the model. Figure 3 
shows the interaction of the GUI with the CAD 
and CLM modules. The main GUI panel 
(Figure 4) manages the input/output information 
from the simulation. The GUI is used to input 
simulation parameters, such as PLA, Mach 
number, and altitude, and output parameters 
from both the controller and CLM. In addition, 
a sub-GUI allows changes to controller, engine, 
and health parameter constants (see Figure 5), as 
well as other parameters that may be varied. 
After running the model, a summary of 
simulation output operating conditions is sent to 
a display in the main GUI. Output that is 
created by the model, along with user defined 
input information is stored in an output structure 
array (a MATLAB data type) that may be 
saved as a MAT-file and reloaded into the GUI 
or accessed independently. 

The components simulated in the MAPSS CLM 
(Figure 6) are the fan, booster, high-pressure 
compressor (HPC), burner, high- and low- 
pressure turbines (HPT and LPT), mixer, 
afterburner (AB), and nozzle. The fan inlet 
temperature and pressure are external inputs 
calculated from user-defined operating 
conditions. High- and low-pressure turbines are 
simulated to extract work from the hot section to 
provide power to the fan/compressor section. 
The CLM also models the forward blocker 
doors, the VABI, the flow through the bypass 



duct, and the bleed flows. The bypass duct 
module determines the pressure, temperature, 
enthalpy, and flow rate up to the point where the 
bypass and core flows mix, known as the mixing 
plane. The CLM does not include an inlet 
model because that is typically not considered 
part of the engine. In addition to these 
parameters, the compressor bleeds and cooling 
flows are simulated. The three state variables 
used in the FORTRAN-based model — low- 
pressure rotor speed, high-pressure rotor speed, 
and the block metal temperature — are also used 
in MAPSS CLM. The overall engine component 
system requires iteration steps to ensure a 
balance of mass flows and energy. 

The CAD module contains the controller and the 
actuator dynamics sub-modules. The controller 
sub-module (Figure 2) uses calculated and 
simulated sensed output values from the CLM to 
determine if the desired set points have been 
reached. The set points are based on the user 
defined operating conditions of PLA, altitude, 
and Mach number. The controller uses maps to 
schedule the actuator demands for the fan, 
booster, and HPC stator positions using the 
engine’s operating conditions. To obtain the 
actuator demands for wf36, a8, al6, the 
controller has five different modes, each of 
which is a proportional-plus-integral (PI) control 
action regulator. The High power mode handles 
CLM operation for PLA > 42.5°. The Medium 
power mode handles CLM operation for 
PLA e (32.5°, 37.5°]. The Low power mode 
handles CLM operation for PLA < 32.5°. As the 
controller transitions between the three power 
modes, percentages of each regulator output, 
based on the current value of PLA, fan stall 
margin, and xnl, are added together to obtain the 
final actuator demand. This is a technique 
known as blending. In addition to the power 
modes there is the Fan Stall Margin mode 
controls the fan in the stall region, and the 
Overspeed mode limits the speed of xnl. These 
last two modes may be turned on and off upon 
request. 

The actuator dynamics sub-module (Figure 7) 
simulates the dynamics of the torque motors, 
servomechanisms and feedback to determine 
actuator positions for the fan, HPC, booster, 


N AS A/TM— 2003-2 1 2543 


guide vanes; nozzle, mixer, and VABI areas; and 
the fuel flow actuator. Lags are used to simulate 
actuator dead-time or phase shift. These lags are 
also captured in MAPSS, even though their 
precise nature has not been explored. 

Implementation Approach and 
Validation 

Before development of this environment 
commenced, a fundamental understanding of the 
FORTRAN engine component models, and the 
information passed between them was obtained. 
This knowledge was used in implementing the 
engine simulation in a block diagram format. A 
script written in MATLAB®, an object-oriented 
software language used for algorithm 
development, modeling, and data analysis, was 
used to take data generated by the FORTAN 
model, and compare it to data generated by the 
block diagrams. The block diagram depicting 
the engine components in MAPSS can be seen 
in Figure 8. This section explains the 
development approach and validation of MAPSS 
then shows example comparisons between it and 
the FORTRAN model. 

As seen in Figure 9, the values that are used to 
determine the fan exit parameters are calculated 
in the intermediate subsystem block diagrams. 
The subsystems contain algebraic equations and 
maps, which are used to characterize the 
behavior of the fan component. For instance, 
the “IGV Scalars” subsystem uses maps to 
determine how the fan inlet guide vane setting 
affects the output of the fan, and the “Corr. 
Airflow w/ Re # Scalars” uses a combination of 
maps and algebraic equations to determine the 
corrected flow and Reynolds number. 
Generally, in developing the subsystems within 
each component model, special attention was 
paid to highlighting the flow of information. 
The general organization of the subsystems in 
this and other component models follows the 
organization of equations in the FORTRAN 
code. However, to simplify the illustration of 
the information flow, the order of execution of 
some of the equations was changed. This would 
not have an affect on the results produced by 
each component. 



The validation process consisted of two steps. 
First, there was an open-loop comparison of 
each block diagram component with the 
FORTRAN model, then a comparison between 
the MAPSS and FORTRAN closed-loop 
models. Each completed block diagram, for 
both intermediate blocks and component blocks, 
was validated before it was linked to another 
subsystem. The MATLAB* script that was used 
for model comparison supplied FORTRAN 
input variables from the subroutines to the block 
diagram and produced comparison plots for each 
output variable. The plots displayed raw data 
from the FORTRAN model and block diagram, 
as well as the percent error between the two, 
which was used to determine the validity of the 
block diagrams. Since the state space equations 
of the FORTRAN routine were exactly 
represented with proper integration, comparison 
utilizing frequency response was not deemed 
necessary and the time-domain error criterion 
was sufficient means of comparison. 

As an example of the comparison plots, Figure 
1 0 shows the response of pnrm2 and w2 for both 
models. The operating conditions for this test 
consisted of a step increase in PLA from 21 to 
30, at sea-level, with a Mach number of 0. From 
these figures, it can be seen that the largest 
difference occurs during the transition from 21 
to 30 (i.e., the transient region). There is very 
little difference (less then a tenth of a percent) in 
both steady state regions. Generally, if the error 
was less than one percent the block diagram was 
considered valid. However, if a block diagram 
had an error that was higher than one percent, it 
was examined to determine the reason for the 
difference. During validation, it was acceptable 
to have an error around one percent during the 
transient, just as long as the steady state error 
before and after was less then one percent. 
According to the plots shown in Figure 10, there 
is good matching of the fan component, 
therefore, it is considered a valid block. These 
steps of block diagram building and validation 
were conducted repeatedly until the CLM and 
CAD modules were fully modeled. 

The CLM and CAD modules were integrated 
(i.e. connected) to form a closed-loop model. In 
order to properly integrate these two modules, 
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the CLM had to be balanced in mass and energy. 
A challenge encountered during implementation 
was developing a scheme that would present a 
balanced model to the CAD module. In the 
FORTRAN simulation, as stated before, a 
while loop is used to ensure the selected mass- 
flow, energy, and work of the turbines is 
balanced out before the next simulation time 
step of 0.02 seconds. For MAPSS, an 
equivalent mechanism would be to feed these 
parameters back into the CLM (i.e. create an 
algebraic loop). However, an algebraic loop can 
prevent the MAPSS CLM from finding a 
convergent solution without the parameters 
being updated. A solution to this problem was 
to construct an iteration block that adjusts a 
vector of parameters, obtained from the CLM, in 
order to drive a norm of component error values 
to zero, and then delay the iteration of the block 
by one simulation time step (40 milliseconds). 
This allowed the components to achieve an 
energy balance within the same tolerances as the 
FORTRAN components, meaning that MAPSS 
CLM and the FORTRAN CLM are running at 
effectively the same sample rate. The adjustable 
vector consists of mass-flow, energy, and work 
parameters that are used to update the CLM at 
each time step, but delayed by one time step to 
avoid an algebraic loop. The error values are 
derived from turbine losses, and booster and 
HPC mass flow rates, as well as bypass, mixer, 
and nozzle throat pressure. They are calculated 
as the difference between the estimated value, 
obtained from component maps, and their 
calculated values. To ensure the balanced 
condition is met before the next control step, the 
iteration block and the CLM, as stated before, 
operate at a higher frequency then the controller. 

Figure 1 1 shows a comparison of state variables 
between the FORTRAN and MAPSS closed- 
loop simulations, for two PLA step up 
conditions. The first is a small step up increase 
from 21 to 22 degrees and the second is a larger 
step up increase, from 21 to 50 degrees. Mach 
number and altitude remained the same as in the 
fan component test. Figures 11(a) and 11(b) 
show that there is very little difference (less than 
one percent) between MAPSS and the 
FORTRAN simulation. The only caveat to this 
result is that, for both models, it was found that 



the bypass ratio at the bypass duct entrance 
(known as bprl4g) became negative during the 
transient (see Figure 12). This result manifests 
itself for PL A >35. This occurs because reverse 
flow is allowed to occur at the bypass duct 
entrance (known as wl4) during the transient. 
Although, through testing this phenomenon does 
not seem to cause any adverse effects in the 
performance of MAPSS, it is undesirable to 
have reverse flow anywhere in the model 
because it may cause instabilities to occur 
among the neighboring components. Therefore, 
the measurements for bprl4g are now prevented 
from becoming negative in the MAPSS CLM 
(as shown in Figure 12), which causes the 
differences between the two simulations shown 
in Figure 13. 

Summary of Results 

This paper presents progress on the development 
of a generic turbofan engine simulation in 
Simulink®. The goal of this effort was to 
develop and demonstrate a flexible simulation 
platform for future research in propulsion 
system control and diagnostic technology. A 
FORTRAN -based model of a modem, high- 
performance, military-type turbofan engine was 
used as the test case for the platform 
development. The focus of this paper was on 
the development of the MAPSS environment 
and its validation using the FORTRAN -based 
engine model. The block diagrams in MAPSS 
were developed from the state space and non- 
linear algebraic equations in the FORTRAN 
model. The validation procedure compared the 
block diagram component models against the 
FORTRAN component models, using a 
MAT LAB script, and then compared the 
MAPSS and FORTRAN closed-loop models. It 
was found through open-loop and closed-loop 
comparison, for a variety of operating 
conditions, that there was less than one percent 
difference between the two simulations when 
looking at the individual components, the 
controller models, and the overall models. 


There was, however, a change made in the 
MAPSS CLM to prevent the bypass ratio from 
becoming negative. Future improvements to 
MAPSS include improvements to the GUI and 
the implementation of a commercial, high- 
bypass turbofan engine in the MAPSS 
environment. 
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APPENDIX A: Figures 



Figure 1: Schematic of turbofan engine model with labeled actuators (above) and sensors (below) [2], 
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Figure 3: MAPSS Module Interaction. 
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Figure 5: MAPSS Change Constants GUI. 
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Figure 8: MAPSS Component Level Model. 
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Figure 9: MAPSS Fan Component Model. 












Comparison for Variable pnrm2 



max(%err)= 0.133879, min(%err)= -0.0317556 



(a) Percent Corrected Fan Rotor Speed 


Comparison for Variable w2 



max(%err)= 0.1 52239, min(%err)= -0.0360609 



(b) Fan Exit Mass Flow Rate 


Figure 10: Open-loop fan exit parameter comparison plots (normalized). 
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Comparison for Variable xnl 



(a) PL A = 21 - 22 deg (Ground Idle Condition) 


Comparison for Variable xnl 



(b) PL A = 21 - 50 deg (Takeoff Condition) 

Figure 1 1: Closed-loop comparison of model state variables (normalized). 
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Comparison for Variable bprldg 



Figure 12: Bypass ratio measurements at forward bypass duct (normalized). 


Comparison for Variable xnl 



Figure 13: Closed-loop comparison of model state variables (bprl4g > 0) (normalized). 
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